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Abstract 

Disturbance of the Maxwellian plasma may occur in the vicinity of a spacecraft due to 
photoemission, interactions between the spacecraft and thermospheric gases, or electron 
emissions from other devices on the spacecraft. Significant non-maxwellian plasma 
distributions may also occur in nature as a mixture of ionospheric and magnetospheric plasmas 
or secondaries produced by photoionization in the thermosphere or auroral precipitation. The 
general formulas for current collection (volt-ampere curves) by planar, cylindrical, and 
spherical Langmuir probes in isotropic and anisotropic non-maxwellian plasmas are examined. 
Examples are given of how one may identify and remove the non-maxwellian components in 
the Langmuir probe current to permit the ionospheric parameters to be determined. Theoretical 
volt-ampere curves presented for typical examples of non-maxwellian distributions include: 
two-temperature plasmas and a thermal plasma with an energetic electron beam. If the non- 
ionospheric electrons are Maxwellian at a temperature distinct from that of the ionosphere 
electrons, the volt-ampere curves can be fitted directly to obtain the temperatures and densities 
of both electron components without resorting to differenting the current. For an arbitrary 
isotropic distribution, the current for retarded particles is shown to be identical for the three 
geometries. For anisotropic distributions, the three probe geometries are not equally suited for 
measuring the ionospheric electron temperature and density or for determining the distribution 
function in the presence of non-maxwellian background electrons. 
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Introduction 

Langmuir probes have been used for many years on rockets and satellites to measure 
ionospheric electron temperature, T e , and ion and electron number densities, /V and N e , see 

references in Boyd 1 and Brace 2 . The Langmuir probe technique involves exposing a metallic 
collector to the plasma surrounding the vehicle, usually by mounting it on a boom that is longer 
than a Debye length at the lowest density expected to be encountered. The boom places the 
probe in contact with ionospheric plasma at a distance where the measurements are undisturbed 
by the presence of the vehicle. The probe current, /, (the sum of the ion and electron currents, 

/, and I e is measured as the collector voltage, V, is repeatedly swept through a suitable range 
with respect to the plasma potential. Analysis of the resulting I-V curves usually yields the 
ionospheric parameters T e and N e . Various sources of measurement error which affect this 
analysis were discussed by Brace 2 , but that work did not cover the effects of non-ionospheric 
or non-maxwellian electron populations. 

The measurements are very simple and direct when the electron energy distribution is 
essentially Maxwellian as is true nearly everywhere in the ionosphere. In this case, the electron 
energy can be characterized by a single scalar value, T e . However, if additional electron 
populations are present in significant magnitude, the determination of T e is more complicated 
and may not be possible if the density of the additional component is sufficiently high. Sources 
of such additional populations may include; (1) photoelectrons produced in the daytime 
ionosphere in the process of forming the ionosphere, (2) degraded secondary electrons 
produced by precipitating auroral electrons and ions, (3) photoelectrons ejected from sunlit 
spacecraft surfaces, (4) secondaries produced by the impact at spacecraft velocities of 
thermospheric molecules, (5) secondaries emitted when the spacecraft potential is high enough 
to accelerate ionospheric ions or electrons to energies that exceed the ionization potential of 
spacecraft metallic surfaces, and (6) electron beams or plumes emitted from devices onboard 
the satellite. The geophysical sources (1 & 2), although present much of the time, usually have 
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a negligible influence on the I-V curves, except perhaps in regions of very low ioilospheric 
density. The non-geophysical components (3-6) can have high enough densities to significantly 
distort the 1-V curves and affect the determination of T t by distorting the electron retarding 
regions of the I-V curves. They may also affect the measurements of total ionospheric density, 
N e and N r In these cases the ionospheric parameters (of the low energy electrons) can only be 
obtained by identifying the non-ionospheric contributions to the I-V curves and subtracting 
them out, or by fitting the I-V curve with theoretical expressions which include the non- 
ionospheric components. Examples of the effects of non-maxwellian situations are: impact 
ionization on the Pioneer Venus Orbiter discussed by Whipple et al. 3 ; an example of two- 
temperature distributions observed in the Venus ionotail described by Brace et al. 4 ; and multi- 
component superthermal electrons by Knudsen and Miller. 5 

This work presents the formulas for and examples of Langmuir probe current collection in such 
non-maxwellian plasmas, and shows when the electron energy distribution can be derived from 
measurements of the first and second derivatives in the retarding region of the I-V curves. 
Specific examples are also presented to illustrate the method of deriving ionospheric parameters 
from fitting the I-V curves. The results show that different probe geometries are affected 
differently by the presence of non-maxwellian components. 

I. The Conventional Langmuir Probe Technique for Maxwellian Electrons 

A Langmuir probe may have any geometry, but cylindrical, spherical and planar probes 
are usually employed because of their symmetry. The basis for T measurements is the 
conventional Langmuir probe theory of Mott-Smith and Langmuir 6 , assuming a Maxwellian 
energy distribution, F(E) , given by: 
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where k is the Boltzmann constant. Thus, in a Maxwellian plasma, when the probe to plasma 
potential, V, is driven negative, / e decreases exponentially 

'. = 4>e* p(%J. (2) 

where / 0 is the random electron current given by A N e e (kT/2n mj' /2 , where e is the 

electron charge, m e is electron mass and A is the probe area. This equation for what is known 
as the electron retarding region in a Maxwellian plasma is the same for all probe geometries. It 
is shown below that this result has a generalization that the expression for the retarded current 
is identical for the three probe geometries for isotropic non-maxwellian distributions. The 
derivation of T e is usually done by fitting the electron retarding region with an exponential 
function, and sometimes T e is obtained electronically by measuring the ratio of the first and 
second derivatives of /, rather than through analysis of the I-V curves themselves. 7 8 

Next consider the ion and electron saturation regions of the I-V curves. The form of the 
acceleration region current depends on the geometry of the collector and is approximated by, 

!+%•)" 0 ) 

where p = 0, 1/2, or 1 for planar, cylindrical, or spherical geometries, respectively. For 
satellite applications Eq. (3) is modified to include the effects of the large spacecraft drift 
velocity relative to the essentially stationary ionospheric plasma. 9 The general ion and electron 
acceleration current expressions are used to generate the I-V curves that are shown in the 
examples rather than the Eq. (3) approximation. 

Examples of the I-V curves for the three geometries in purely Maxwellian plasmas are given in 
Figures la-c which show I-V curves for 4 different temperatures. The calculation assumed N e 
= 10’’ cm ’, a collector area of 6.73 cnr and a spacecraft velocity of 8 km/s. Note that the curves 
are down-ranged by a factor of 10 where necessary to allow the large electron acceleration 
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currents to be shown on the same plot us the rest of curve. These figures illustrate'the effect of 
T, on the width of the electron retardation regions. 

For a Maxwellian plasma, all three probe geometries are suitable for the measurement of T e . 
However, the electron saturation regions differ greatly with collector geometry. Note that the 
electron saturation region of the cylindrical probe is nearly independent of T e , whereas this 
region is highly temperature-dependent for the planar and spherical probes. This illustrates a 
practical advantage of the cylindrical probe; i.e., that a fixed probe potential can be applied to 
make continuous measurements of N e without pausing to measure T e , assuming a known and 
stable value for the potential of the spacecraft. The cylindrical and spherical probes have the 
advantage of producing large saturation currents for the same density, an advantage for 
measurements in regions of very low density. 

II. Multi-component Plasmas and Other Non-maxwellian Distributions 

When non-ionospheric (non-maxwellian) electrons are present in detectable quantities the 
standard Langmuir equations may lead to significant errors in the measurements of N e and T e . 
To handle these situations we examine the current equations for arbitrary electron energy 
distributions for the three probe geometries. These equations are integrals over the velocity 
distribution functions (VDF) for multi-component populations of ions and electrons. The 
discussion begins with formulas for general anisotropic VDFs, then considers the general 
isotropic case, and finally the anisotropic case for a Maxwell-Boltzmann distribution. The 
effects of these non-ionospheric populations are illustrated by calculating I-V curves for a two 
temperature distribution as an example of the isotropic case, and an energetic electron beam 
with a superposed Maxwell-Boltzmann distribution as an example of the anisotropic case. The 
information that can be obtained about the electron distribution functions by using the 1st and 
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2nd derivatives of the I-V curves is discussed, and later fits are made to the curves to illustrate 
obtaining the ionosphere parameters. 

III. Current Equations for Arbitrary Velocity Distributions 


The general form of the VDF for the sth specie (s is an electron or an ion specie) divided 
into distinct populations j is a sum over the j populations: 

F s («, v,w) = £ n i F s ( M - v ’ w ) ’ ( 4 ) 

j 

where n J s is the number density of population j for specie s, Fj is the individual VDF, and 
u,v,w are velocities in the x,y,z directions. The sum of the number densities over the j 

populations is the total number density of specie s, ^^n J s = n s , and the VDFs are normalized 

j 

to unity. 


oo oo oo 

u, v, vv) = 1 . 

— oo — oo — oo 



(5) 


Two classes of distributions: isotropic, and anisotropic are considered. For the isotropic VDF, 
the general form is, 

F s ( w - V, w) isolropic = Gi (u 2 + V 2 + w 2 ) , (6) 

where G is a function of the sum of velocities squared or of the energy. The most applicable 
form for a collision dominated plasma is the Maxwell distribution function: 



where m and T are the mass and temperature of the jth component of the sth specie (for 
simplicity, the subscripts and superscripts are not explicitly written in all cases). 
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The anisotropic distribution has a preferred direction in space and could result, for 
example, from a beam of electrons induced by an acceleration process, or from rapid motion of 
the probe through the plasma, or from an anisotropy in the temperature of the electrons or ions 
(electrons parallel and perpendicular to the magnetic field can have distinct temperatures). 
Anisotropic T e was first measured by Clark et al. lu and theoretically predicted by Schunk and 
Watkins." We will employ a specific form of the anisotropic VDF which applies to many 
space plasma situations: 

F s (“’ V ’ W ) anisotropic = F s («(“ ~ «0 )A V ~ V 0 )’ C { W ~ "0 )) isotropic > ( 8 ) 

where a,b, and c are constants that could represent temperature anisotropy, and u 0 , v 0 , vv 0 

represent a drift velocity of the plasma relative to the probe. An applicable form of the 
anisotropic VDF is the Maxwell-Boltzmann distribution which is given by Eq. (7) with the sum 
of squares of the velocities replaced by the sum of squares with respect to the drift velocity, 

a 2 (u-u 0 f + b 2 (v-v 0 )~ +c 2 (w- w 0 ) 2 . 

Another useful form for the anisotropic VDF is as a spherical harmonic expansion. 

III.A. General Current equations 

The Langmuir probe current equations are given below for the three standard geometries; 
planar, cylindrical, and spherical. To simplify the formulas we assume large symmetrical 
sheaths and orbital-motion-limited collection for the cylindrical and spherical probes, which is 
valid over essentially all ionosphere conditions as long as the relevant probe dimension is 
smaller than a few cm. It will become apparent that each geometry has advantages and 
disadvantages which depend upon the type of non-maxwellian distribution. 

III.A. 1. Planar Probe General Formulas 
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For a planar probe of area A, the general formula for the current is. 


oo 



where q is the charge, in the mass of the sth specie and jth population. The lower limits on the 
integral correspond to accelerated (qV> 0), retarded (c/V<0) particles respectively, it is the 
component of velocity perpendicular to the planar probe. And Fid is the one dimensional or 
"ID" distribution, defined by the integral over velocities parallel to the probe surface. 


F{u) w = | jdvdwF(u,v, vv), (10) 


and where the subscripts and superscripts are dropped for simplicity of notation. This result is 
general for any VDF, isotropic or anisotropic. In the ionosphere, the electron current at zero 
probe potential is about a factor of 170 greater than the ion current because of the large 
ion/electron mass ratio, hence the I-V characteristic is dominated much more by the electron 
current than the ion current. Therefore, the total current in the retarding region can be used to 
attempt to determine the electron velocity distribution. We adopt the convention that the electron 
current is positive. Thus the 1st derivative of the planar electron current is given by, 

_ Ae Y n >F' 
deV m e j 

where eV<0. The derivative is proportional to the "ID" VDF. This result does not imply that 
the individual populations) can be obtained from the 1st derivative which gives the sum of the 
populations. However, if the populations dominate in distinct energy regions, then the 
derivative can be used to determine the dominant populations. Also the “ID” VDF is not the 
same as the original “3D” VDF of Eq. (4), however, it is a useful concept. Later in the 
discussion on isotropic plasmas, the relationship between the “ID” and “3D” distributions is 
examined. In general, the 2nd derivative of the planar electron current appears to be not related 
directly to the distribution function, since it is the 1st derivative of the "ID" VDF. However, if 


-2eV 


in 


(ID 
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the distribution takes the specific form of anisotropic VDF that we assumed in Eq. (8), then it 
can be shown using integration by parts, that. 



where a,b,c are anisotropy parameters which are distinct for each population, and u o is the 

population dependent drift velocity in the direction perpendicular to the planar probe surface. 
Thus the 2nd derivative of the planar current is equal to a weighted sum of the "3D" VDFs for 
our assumed form of the anisotropic distribution. If the populations have distinct energy 
domains, then it may be possible to obtain some information on the VDF of each of the 
populations. The result is independent of the drift velocities in the directions along the probe 
surface, consistent with the planar probe sensing the ”1D” distribution. The planar probe 
current was treated by Federov 12 for the anisotropic distribution expanded in spherical 
harmonics. In the special case of an isotropic VDF, (a=b=c= 1, u= 0) then Eq. (12) reduces to 
the sum over the individual "3D" VDFs. 

III.A.2. Cylindrical Probe General Formulas 

For a cylindrical probe with an anisotropic VDF, the general formula for the current is 
more complicated than for a planar probe. The geometry is no longer rectilinear and we must 
integrate over the circumference of the cylindrical surface which changes the direction of the 
normal and tangential velocities relative to the anisotropy direction in space. If we let the z-axis 
be along the probe axis, the dependence on a drift velocity, w Q in the z direction is eliminated. 

We use our general form for the anisotropic VDF, given by Eqs. (6) and (8) and define a "2D" 
VDF, 
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F(a(u-u {) ),b(v-v 0 )) if} = jdwF(u,v,w) , (13) 

itnisotnipic 

where F tmixilln)pit . is given by Eq. (8). At an arbitrary position 0 along the probe circumference, 

the radial, u , and tangential, v , velocities are related to the x,y components u.v by the 
transformation, 

u = u cosQ - v sinQ, v = u sin 9 + v cosQ. (14) 

The orbital-motion-limited cylinder probe current is then given by, 

In °° 

icy Under = 1 ^ J n Ws 2n j dd J ud “j“ 2 + Fj (a(u cos 9 - u 0 ),b(usin 9 - v 0 )) 2D , (15) 

S ' J 0 o, 

V m 

where A is the cylinder area and transformation ( 14) was used in Eq. (13), u is relabeled u and 

we may set v =0, Mott-Smith and Langmuir. 6 The integral is complicated and involves 
integration over a modified Bessel function in the case of a Maxwell-Boltzmann VDF. The 1st 
and 2nd derivatives of the cylindrical probe current in general do not directly give the 
distribution function. However, Federev 12 has shown that the cylindrical I-V curve can yield 
information on the anisotropic distribution function. 

III.A.3. Spherical Probe General Formulas 

For the spherical Langmuir probe, there is no preferred direction and the current is an 
integral over 471 directions on the sphere. The general transformation of velocities from the 

fixed x,y,z velocities, u',v',w' to velocities on the sphere at the position given by the spherical 
polar angles, 9, 0 with u the radial velocity and v and w the tangential velocities, is given by, 

u = u cos 9 + sin cos 0 + w sin 0) 

v = —u sin 9 + cos 9(v cos 0 + vv sin 0) (16) 
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w = - v sin 0 + vv cos 0 

where 6, 0 are the polar and azimuthal angles respectively. The inverse transformation is the 
transpose ot (18). Then the spherical probe current in the limit of a large sheath (orbital- 
motion-limited case) is given by, 5 

n In °o 

‘sphere = J ^ n s4s J sin 6cW J J U(iu[u 2 + • 


s-J 


0 -2 qV 


f M w cos 0 - m 0 ),/?(« sin 0 cos 0- v o ),c(wsin0sin0 - w 0 )) , (17) 

where A is the sphere area, and F is our general anisotropic "3D" VDF given by (8). For 
retarded electrons, the 1st derivative is, 


di. 


sphere 


deV 


a it In <» 

= — ^n J e e-^ jsin Odd j dip | uduF 1 , 


(18) 


2fV 




and the 2nd derivative is given by, 


d 2 i 


deV 




[^277 - 


(19) 


Thus the 2nd derivative of the spherical probe current is the average over all orientations of the 
anisotropic distribution function, or it yields the isotropic part of the distribution as was shown 
by Federov 12 . 


III.B. Current Equations for Isotropic Distribution of Electrons 

The isotropic velocity distribution function has the general form given by Eq. (6). We now 
demonstrate that, for isotropic electron populations, the three geometries have I-V curves with 
identical functional form for the retarded current. Thus for all three geometries, the 1st 
derivative yields the “ID” VDF and the 2nd derivative yields the “3D” VDF. This result means 
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that theoretically all three are suitable for measuring isotropic, non-maxwellian electron velocity 
distributions. 

III.B.l. Standard Form of the Isotropic, Retarded Current 

Given the form of the isotropic VDF, Eq. (6), the “ID” distribution can be rewritten as, 

oo 

G w (u 2 ) = 2KjtdtG(u 2 +t 2 ), (20) 

o 

where G is the “3D” VDF and the polar transformation of variables was used, v = t cos a , 
w = rsina . Thus the planar current can be rewritten, 

oo oo 

ip° Irop,c = A^^nq J udu27t ^ tdtG{u 2 + t 2 ^j , (21) 

o.V^ o 

O T/ 

where X = ; the lower limit 0 is for accelerated particles when X>0, and the limit V — X is 

m 

for retarded particles when X<0. The accelerated current is independent of voltage, and thus we 
consider only the retarded particles which yield information on the VDF. With a change of 
variables, the retarded planar current expression is, 


jiso tropic 
P retarded 




( 22 ) 


The cylindrical current is obtained by substituting an isotropic VDF in Eq. (15), giving the 
form. 
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OO OO 

■noliopn _ 2iuj J* ttdu'Jir + X J dwG(ll~ + H’ 2 j . 


(23) 


(W-X 


Successive transformation of variables, u 2 + X = r, and t = sin0j\ w = cos ds, and the 
integral over 9 yields 4 , and thus the retarded cylindrical current is given by the integral in 


Eq. (22), 


etarded p 


__ • isotropic 


retarded * 


The spherical current, Eq. (17), for an isotropic distribution has the form, 


^tropic = A ^nqK J udu(u 2 + x)g(u 2 ), (24) 

o.d^x 

and using the transformation, u 2 + X = s 2 , the spherical probe retarded current is found to 
have the form of Eq. (22), or, i‘ s sotropic ' retarded = i l *° ,ropic * ^ ^ . Thus the three probe 
geometires have the identical retarded current expression for an isotropic VDF. 


III.B.2. First Derivative of the Isotropic Retarded Current 

The 1st derivative of the retarded current was shown in Eq. (1 1) to be given by the “ID” 
distribution for a planar probe, thus for isotropic plasma, the three probe goemetries all have 
the identical form given by, 


j: isotropic 

deV 




( 2eV ' 


'ID 


m 


(25) 


e J 


where the “ID” VDF is given by Eq. (20). 

Examine the relation of the “ID” distribution to the full “3D” VDF. The “3D” VDF can be 
expanded in Laguerre polynomials 13 , 
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oo > 

c(. r)=<"J]nfW, (26) 

n=() 

then the “ID” VDF is the integral, 

oo 

G iD {x) = /rj dsG(s + jc), (27) 

o 

which can be rewritten, using the relations among Laguerre polynomials given by Rainville 13 , 
as follows, 

oo 

G w (x) = e-*KY,gA'\x)- ( 28 ) 

n = 0 

For a Maxwell distribution g n =0 for n>0, which is the only case when the “ID” and “3D” 
distributions have the same functional form. The first few terms of the two distributions are 
given by, 

G{x) = e~ x S\g 0 +g l (l-x) + g 2 (l-2x + x2 ^ + ..)j (29) 

G w {x) = k <T- ( jgo - g\x + g 2 [-x + x / 2 ^j + ...j. (30) 

The distributions are distinct, but the dominant terms in jc are identical. 

III.B.3. Second Derivative of Isotropic Retarded Current 


The 2nd derivative of the isotropic current is readily obtained from Eq. (12) by 
substituting the isotropic distribution, we obtain, 


d l [isotropic 2k 


deV- 


n 0 eG\ 


f -2eV^ 


n 


V m e ) 


( 31 ) 


A result first found earlier by Mott-Smith and Langmuir. 6 In the next section we give the 
current formulas for the special anisotropic case of the Maxwell-Boltzmann distribution with 
isotropic temperatures. 
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III.C. The Maxwell-Boltzmann Anisotropic Case 


The VDF for a Maxwell-Boltzmann (MB) distribution is given by the formula, 


F s ! (u, v, w) MB = 


mi 


2nkT< , 


exp 


mi 


2kTj 


J ((« - «o ) 2 + f + i w ~ w 0 f ) 


(32) 


where the drift velocity components u fJ v w g may be different for each specie and population 
of a specie. For simplicity, we will assume below that all species have the same drift velocity. 


III.C. 1, Planar Probe Equations for Anisotropic MB Distribution 


The "ID" VDF is obtained by substituting (34) in Eq. (12), 

u-u Q f , 

) 

which leads to the current expressions, 

1 planar MB accelerated = X ' random ( eXp (~ r ) + V?T rerfc{-r)} 

j.s 

for accelerated particles, and 

‘P lallar MB retarded = ‘ random ^ eX P^"~( - r ) ) + ^ ^// C (V~*7 ~ r )j * 

j.s 

for retarded particles, where r is the ratio of drift velocity to thermal velocity, 



and 77 is the ratio of voltage to thermal energy, 


fM 


mi 


MB\ 


ID 


,2 mi j 


\Vi r 

exp 


V 


2 kTj 


(33) 


(34a) 


(34b) 


( 35 ) 
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ni-aZ. 

IcTj 


(36) 


and the random current, the current due only to thermal motion of the particles, is. 


, _ . j , AT/ 

1 randoms 


\2nm ] s 


(37) 


III.C.2. Cylindrical Probe Equations for Anisotropic MB Distribution 


Using the Maxwell-Boltzmann distribution (32) in the expression for the cylindrical probe 
current, (15), and transforming the drift velocities, we find, 


( cylinder = ^ ‘random ^ J tdt^t 2 + T] exp(-(r 2 + r 2 ))/ 0 (2/r), 

sJ O.^rj 


(38) 


which formula cannot be evaluated in simple analytical form. Io is the modified Bessel function 
and r and 77 are given above. Useful asymptotic approximations have been obtained for this 
integral by Hoegy and Wharton. 9 However, since we are evaluating the cylinder current for a 
wide range of values of r and 7) , we evaluate it numerically using the fast, efficient routines 
from Numerical Recipes. 14 The limit of Eq. (38) for zero drift velocity or the isotropic Maxwell 
distribution gives for accelerated particles, 

2 


1 random 


-j=^ri + exp(rj)erfc(4ri) 

for large 77 ,this expression is approximated by ■ 

For retarded particles the current is given by, 


(39a) 


random 


exp(rj). 


(39b) 


III.C.3. Spherical Probe Equations for Anisotropic MB Distribution 
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The Maxwell-Boltzmann distribution (32) in the expression for the spherical probe 
current, (17), gives the expression, 

oo 

1 sphere = ^ 1 random ^ J dt{t 2 + o)(exp(-(f - rf ) - exp(-(f + r) 2 ]]. 

S 'J o.4~n 

The spherical geometry, results in spherical Bessel functions which can be integrated to give, 

'sphere accelerated random 

s,j 

for accelerated particles, and 

^P he re re , a rded ~ ^ ‘random ^\( 7 ? + r + + r ) ~ er f(4~ r l ~ r )} + 

s,j 

+ - rf ) - — ^ — exp^-^-sJ-rj + rf 

for retarded particles. In the limit of zero drift velocity, or an isotropic Maxwell distribution, 
Eq. (40a) reduces to, 

uJ'i+i)- ( 41a > 

and Eq. (40b) reduces to, 

e-w( ri)- (4 lb) 

To aid the reader in visualizing the effects of non-maxwellian electron distributions on the I-V 
curves. Sections IV and V show theoretical I-V curves for mixtures of electron populations 
having different temperatures and density ratios. I-V curves are also shown for mixtures of 
thermal electrons and energetic beam electrons. 

IV. I-V Characteristics for Two-Temperature Distributions 



rexp(-r 2 )- 




erf{r) 


(40a) 


17 



10/23/98 3:42 PM 


Figures 2 a, b, and c show planar, cylindrical and spherical probe I-V curves' for various 
density ratios of 10,000 K and 1000 K electron populations, with a combined total density of 
10 5 cm ’. Curves are shown for energetic component densities relative to the total density of 0, 
20, 33.3, and 50%. 0% corresponds to the complete absence of the higher temperature 
component. 

The higher temperature component enhances the current at greater retarding voltages as 
seen by comparing these figures with Figures 1 a-c. The temperatures and densities of the two 
components can be obtained by fitting a two-temperature exponential to the curves or by fitting 
the near and far retarding regions with single-temperature exponential. The latter procedure was 
used successfully in the fitting of I-V curves from the PVO Langmuir probe which exhibited 
two-temperature characteristics in certain regions of the Venus ionotail. 4 

V. I-V Characteristics for Thermal Electrons Plus an Energetic Beam 

Figures 3 a-c show theoretical I-V curves for combinations of 1000 K thermal electrons 
and a 1 eV energetic beam having a 2000 K thermal spread. The assumed total density is 10 5 
cm' 3 . The beam is directed normal to the planar probe surface and perpendicular to the 
cylindrical probe axis. The beam direction is irrelevant for the spherical probe. This case is our 
example of an anisotropic distribution. The non-exponential nature of the curves shows 
immediately that the energy distribution is non-maxwellian. 

These examples are illustrative of the effects of non-maxwellian plasmas on Langmuir 
probe measurements. Other combinations of secondaries, energetic beams and ionospheric 
electrons would change the shapes and amplitudes of the Langmuir probe characteristics in 
different and unique ways. The question remains as to how the ambient ionospheric 
temperature and density can be derived when such additional electron populations are present 


18 



10/23/98 3:42 PM 


as an alternative to the derivative technique. The next section illustrates what can b'e learned 
about such complex distributions by fitting the I-V curves using either one or two temperature 
distributions. 

VI. Fitting I-V curves for one or two temperature Maxwellian distributions. 

If one knew the form of the energy distributions of the secondary electron populations, 
one could in principle fit the I-V curves for the energy and density of each of the contributions. 
In practice this is difficult, since one usually only has the I-V curves themselves to work with. 
After finding that a single temperature fit fails to achieve a sufficiently small standard deviation, 
one could proceed by trial and error to add various other components to the fit to improve the 
quality of the fit. The first step is usually to add another maxwellian component at a suitable 
temperature and density. If such a two-temperature distribution fails to achieve an acceptable 
fit, one can proceed to add or substitute directed energetic beams with appropriate energies and 
temperatures of the type illustrated above in Sections V. 

In this section we go part way down this path by showing how well single and two- 
temperature probe theories fit the I-V curves shown in Figures 1 - 3. We fitted the theoretical 
curves shown in Figures 1-3 using techniques similar to those used in the analysis of Langmuir 
probe measurements from Pioneer Venus Orbiter 4 and Dynamics Explorer-2. 15 The retarding 
regions were fit first to a single Maxwellian. In some cases, the curves were then fit using a 
two-temperature Maxwellian to explore the feasibility of obtaining information on the nature of 
the distributions and to recover the ambient parameters. 

The single temperature fits were made using the simple form: 
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/ = a + bV + cexp 


1 1604.45 V ^ 

71 J 


(42) 


where V is the probe voltage in volts, T, is the single temperature in K, and a, b and c are 
constants to be determined by the fit. The first two terms represent the ion current amplitude 
and slope, respectively. The third term represents the retarded electron current. The fitting 
procedure steps through a range of temperatures T, while performing a least square fit of Eq. 
(44) to the theoretical curves to determine the coefficients. A Newtonian scheme is used to 
refine the fit to obtain a temperature that minimizes the standard deviation, STD. 


A similar method was used to fit the curves for two temperatures using the formula: 


/ = a + bV + c exp 


1 1604.45 V 


3 


+ d exp 


^1 1604.45 


(43) 


The ratio of coefficients d and c determine the ratio of the densities of the high and low 
temperature components. The domains of the two temperatures are kept distinct to ensure the 
uniqueness of the fit parameters. In a flight application, where perhaps millions of curves may 
have to be fitted, a faster scheme would be devised. The present fit procedure, however, is 
adequate to demonstrate the nature of the two-temperature fitting procedure. 


These fit procedures were applied to the three cases of different probe geometry, planar, 
cylindrical, and spherical, and to three distribution functions, purely Maxwellian, two- 
temperature Maxwellian and a Maxwellian with a superimposed electron beam. When the two- 
temperature fits were applied to the one-temperature curves, the density derived for the second 
component was essentially zero, thus confirming that only a single temperature was present. 
The fits to the single Maxwellian are not shown in the table because they reproduce the input 
values of T t and N e with greater than 1 % accuracy corresponding to STDs of less than 5x 10 s . 
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Table 1 lists the derived temperatures and fit STDs for a two-temperature distribution and 
tor ditterent density ratios of the two components. The numbers in parenthesis are power of 
10, e.g., 3.2(-5) = 3. 2x10' \ The curve to be fitted were obtained using the appropriate 
equations that were presented earlier. We assumed no digitization errors and used 226 voltage 
values between -3 V and 0 V. The theoretical temperatures were 1000 K and 10,000 K, 
respectively. Although not shown in Table 1, the two-temperature fits retreived the two 
temperatures and three density ratios with negligible error. The ratio of fit parameters, d/(c+d) 
for the two temperature least square fit were, 0.44, 0.61, and 0.76 which agree with the 
theoretical values corresponding to density ratios njn, of 0.25, 0.5 and 1.0 respectively. These 
ratios are related by the formula. 


d / 

d + c 1 + /’ 


where 



Table 1. T e and STD for single and two-temperature fits to curves produced by a two- 
temperature distribution at 1000 K and 10,000 K. 


ratio n,/n, 

0.25 

0.5 

Planar Probe 

1.0 

Te (single T fit) 

2130 

3694 

5652 

STD (single T fit) 

5(-2) 

7- (-2) 

5. (-2) 

STD (two T fit) 

8.(-6) 

7.7(-6) 
Cylindrical Probe 

9.3(-5) 

Te (single T fit) 

2133 

3698 

5655 

STD (single T fit) 

5.4(-2) 

7.2(-2) 

5.4(-2) 

STD (two T fit) 

8.(-6) 

7.7(-6) 
Spherical Probe 

7.2(-6) 

Te (single T fit) 

2117 

3678 

5633 

STD (single T fit) 

5.4(-2) 

7.3(-2) 

5.5(-2) 

STD (two T fit) 

3.6(-5) 

3.6(-5) 

3.8(-5) 


Note that the STDs of the single temperature fits are much higher than those obtained 
when the more appropriate two-temperature fits are performed. All three geometries provide 
very low STDs for the two temperature fits, indicating that all are appropriate for measuring 
two-temperature electron populations. Table 1 also demonstrates the vast improvement in the 
STD when a two-temperature fit is used instead of the simple single temperature fit. We find 
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that a two-temperature fit still works when the original distribution contains only one 
temperature component. The STD does not improve significantly relative to a single 
temperature fit. This result suggests that the presence of a two-temperature plasma can be 
confirmed by comparing the STDs for a two-temperature fit and a single temperature fit. A 
single temperature plasma is present if the STDs are low for both types of fit. A two- 
temperature plasma is present if the single temperature fit yields a high STD and a two- 
temperature fit yields two temperatures, both with small STDs. If neither of these outcomes 
occur, then a more complicated energy distribution must be present. 


Next it is shown what happens when one of these more complicated energy distributions 
is present. Table 2 shows the T e and STD results obtained when single and two-temperature 
fits are performed on the curves shown in Figures 3 a-c. These curves were derived for the 
case where both a directed electron beam and a single temperature plasma are present. The 
thermal plasma has a temperature of 1000 K, while the beam was given a directed energy of 1 
eV and a superposed temperature of 2000 K. The family of density ratios are as shown in the 
figures. 


Since the curves do not appear as simple exponentials, it is not surprising that the STDs 
are high for both single and two-temperature fits. The fit data for the planar probe are not 
shown because the STDs are so large as to make the results meaningless. 


Table 2. T e and STD from single and two-temperature fits curve produced by a 
Maxwellian component at 1000 K and a 1 eV beam with a superposed 2000 K temperature 
n 2 /n, 0.25 0.5 1.0 

Cylindrical Probe 

T/STD (single T fit) 4935/0.14 1 1,426/9.8(-2) 22,260/5. (-2) 

Two temperature fit yields T,= 1000 K, T,=36,800 K 

STD (two T fit) 4.8(-3) 1.3(-2) " 3.(-2) 

d/(c+d) (two T fit) 0.865 0.94 0.98 

Spherical Probe 

T/STD (single T fit) 3100/8.8(-2) 5726/8.(-2) 8483/4.(-2) 

Two temperature fit yields T^IOOO K, T,= 12,000 K 

STD (two T fit) l.l(-3) 3. (-3) ' 6.7(-3) 

d/(c+d) (two T fit) 0.57 0.75 0.87 
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Discussion 

While this work has presented and illustrated the Langmuir probe theory for non- 
maxwellian electron populations, we have also tried to convey some sense of how one might 
proceed to use these results in the analysis of experimental I-V curves obtained in space. In our 
experience, the vast majority of such I-V curves can be fitted very well using the original 
Langmuir probe theory for a single temperature and density plasma. Under certain conditions, 
however, the effects of photoelectrons or impact secondaries can be seen in the I-V curves as a 
high energy tail in the retarding regions. 4 Although we have not experienced it in the 
ionosphere, one can be sure that a Langmuir probe would be sensitive to the presence of an 
energetic electron beam, and/or any secondaries that it might produce. 

The presence of secondaries is not usually a surprise to the experimenter who understands how 
spacecraft interact with the atmosphere and knows what other devices onboard may produce 
interferring plasma populations. Therefore, one starts off with a good idea of which kinds of 
distributions to employ in fitting such curves. In our experience most such curves can be fitted 
by introducing into the fit a second Maxwellian component. If the spacecraft is highly positive, 
or if an electron beam source is known to be onboard, it may be necessary to add or substitute 
an electron beam term into the fitting routine. In principle, some such combination of 
Maxwellian and beam sources can reproduce almost any I-V curve that is likely to be 
encountered. However, the parameters derived from the I-V curves may not always give either 
the original ionosphere component or the secondary components. Interactions among the 
various populations could conceivably produce a final distribution at the probe location that is 
not a simple superposition of the components that produced it. A detailed evaluation of this 
possibility would be required to validate that assumption in each case. 
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We showed that all three probe geometries are not equally suitable for the measurement of 
T„. An earlier work by Hoegy lh showed how different measurement techniques, for example, 
radar backscatter, a.c. mode Langmuir plate, etc. using distinct methods of deriving temperature 
tend to give different temperature values when the distribution is non-maxwellian. Thus distinct 
techniques or instrument type also lead to different plasma probe results. The figures 
demonstrated that electron saturation regions differ greatly with collector geometry. The 
saturation region of the cylindrical probe is nearly independent of T e , whereas this region is 
highly temperature-dependent for the planar and spherical probes. We showed that the 
retarding region of the I-V curves is identical for all three probe geometries when the 
distribution is isotropic. The 1st derivative of the retarding regions gives the "ID" VDF while 
the 2nd derivative yields the full "3D" VDF. When the energy distribution is anisotropic, it is 
not always possible to derive the distribution uniquely, but some information can be obtained 
from the derivatives. For example, the 1st derivative of the planar probe current and the 2nd 
derivative of the spherical probe current give averages of the anisotropic distribution functions. 
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Figure Captions 

Figure 1 .a. Planar probe curves computed for a Maxwellian plasma illustrate the effect of T e on 
the width of the retarding region. Note that the electron saturation region is flat because the 
calculation assumes an infinite plane; i.e., a perfectly guarded planar probe. The amplitude of the 
electron saturation current varies greatly with T t 

Figure l.b. Same as Fig. I, a. but for a cylindrical probe. Note that the electron saturation 
currents are nearly independent of T e . 

Figure l.c. Same as Figure l.a. but for a spherical probe. 
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Figure 2. a. Planar probe curves for mixtures of electron populations having temperatures 
of 1000 K and 10,000 K, The curves are labeled with the percentage of the total density 
represented by the higher temperature component. The low temperature component 
dominates the shape of the curves at low retarding potentials, while the high temperature 
component dominates the more negative end of the retarding region. 

Figure 2.b. Same as Figure 2a, except for a cylindrical probe. 

Figure 2.c. Same as Figure 2a, except for a spherical probe. 

Figure 3. a. Planar probe curves for a 1 eV energetic beam with a 2000 K energy spread 
and a 1000 K thermal component. The curves are labeled with the percentage of the total 
density represented by the beam. Note that the beam component dominates the retarding 
region for relatively small percentages of beam density, thus potentially interfering with the 
measurement of the cold ionospheric component. 

Figure 3.b. Same as Figure 3 a, except for a cylindrical probe. The retarding region is not 
as sensitive to the beam component as was the case for the planar probe (Fig. 3a). Also 
note that the electron saturation currents to a cylinder are relatively less sensitive to the 
presence of a beam than are planar probes and spherical probes (Fig. 3c). 

Figure 3.c. Same as Figure 3a, except for a spherical probe. The retarding region is less 

sensitive to the beam component than either the planar or cylindrical probe (Figures 3a and 3b) 

because thermal electrons are collected from 471 radians, while the beam approaches the 
collector from only one direction. 
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Planar Probe, Maxwellian Distribution 



Figure 1 .a. Planar probe curves computed for a Maxwellian plasma illustrate the effect of T e 
on the width of the retarding region. Note that the electron saturation region is flat because the 
calculation assumes an infinite plane; i.e., a perfectly guarded planar probe. The amplitude of the 
electron saturation current varies greatly with T e 


Cylindrical Probe, Maxwellian Distribution 



Figure l.b. Same as Fig. l.a. but for a cylindrical probe. Note that the electron saturation 
currents are nearly independent of T e . 





Spherical Probe, Maxwellian Distribution 



Figure l.c. Same as Figure l.a. but for a spherical probe. 


Planar Probe, Two Temperature Distribution 



Figure 2.a. Planar probe curves for mixtures of electron populations having 
temperatures of 1000 K and 10,000 K. The curves are labeled with the percentage of the 
total density represented by the higher temperature component. The low temperature 
component dominates the shape of the curves at low retarding potentials, while the high 
temperature component dominates the more negative end of the retarding region. 





Cylindrical Probe, Two Temperature Distribution 



Figure 2.b. Same as Figure 2a, except for a cylindrical probe. 


Spherical Probe, Two Temperature Distribution 



Figure 2.c. Same as Figure 2a, except for a spherical probe. 





Planar Probe, Energetic Beam Distribution 



Figure 3. a. Planar probe curves for a 1 eV energetic beam with a 2000 K energy spread 
and a 1000 K thermal component. The curves are labeled with the percentage of the total 
density represented by the beam. Note that the beam component dominates the retarding 
region for relatively small percentages of beam density, thus potentially interfering with the 
measurement of the cold ionospheric component. 


Cylindrical Probe, Energetic Beam Distribution 
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Figure 3.b. Same as Figure 3 a, except for a cylindrical probe. The retarding region is 
not as sensitive to the beam component as was the case for the planar probe (Fig. 3a). Also 
note that the electron saturation currents to a cylinder are relatively less sensitive to the 
presence of a beam than are planar probes and spherical probes (Fig. 3c). 





Spherical Probe, Energetic Beam Distribution 



Figure 3.c. Same as Figure 3a, except for a spherical probe. The retarding region is less 
sensitive to the beam component than either the planar or cylindrical probe (Figures 3a and 
3b) because thermal electrons are collected from 4 n radians, while the beam approaches the 
collector from only one direction. 




